!===============================================================
!
! FUNZIONE SCHUBERT
!
!===============================================================
!			VARIABLE BOUNDS
!===============================================================
!       -1.D+1               1.D+1         Limits of variable 1
!       -1.D+1               1.D+1            "        "      2

SUBROUTINE SCHUBERT(X,N,F)

      IMPLICIT NONE

      INTEGER          :: N
      DOUBLE PRECISION :: X(N), F


      DOUBLE PRECISION :: F1,F2
      INTEGER          :: I

      F1 = 0.D0
      F2 = 0.D0

      DO  I=1,5
        F1 = F1 + DBLE(I)*DCOS((DBLE(I)+1.D0)*X(1)+DBLE(I))
        F2 = F2 + DBLE(I)*DCOS((DBLE(I)+1.D0)*X(2)+DBLE(I))
      END DO

      F = F1*F2

	  !F = F + 186.730908831024d0

      RETURN

END
